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VHF radar echoes are greatly enhanced near the zenith relative to other 
directions. This enhancement must be due to reflection from horizontally 
stratified laminae of refractive index. In this paper, we suggest that the 
refractivity laminae are due to the displacements of low frequency buoyancy 
(internal gravity) waves acting on the background vertical gradient of 
refractivity. 

The radar cross section <j(k) is given by (OTTERSTEN, 1969) 
a(k) = j k 4 $u(k) 

where k •» 4ir/X adar > V 1 is the potential radio refractive index, and 0p()^) is 
the spatial powlrspectral density of u. If_the fluctuations of p are due to a 
spectrum of vertical displacement acting on M = <dp/dz>, the background 
gradient of p, then 

$p(k) = M 2 K r (k) 

«\, £ 'Xj 

where Eg (Jc) is the spatial power spectrum of vertical displacement. In order 
for o(kJ to be strongly enhanced near the zenith, E (Jc) must also be strongly 
enhanced. 


VANZANDT (1982) has shown that the observed spectra of mesoscale wind 
fluctuations in the troposphere and lower stratosphere can be modeled by a 
universal spectrum of buoyancy (internal gravity) waves that is a slight 
modification of the GARRETT and MUNK (1975) model of oceanic internal gravity 
waves. Since the observed frequency spectrum is red, the buoyancy wave model of 
the vertical displacement spectrum is strongly enhanced near the zenith. In 
other terms, the resulting refractivity irregularities are strongly stratified. 
The model spectrum is 


E (k) 
C 



where n is the buoyancy (Brunt-Vaisala) frequency, E^ is the normalized energy 
per unit mass (E^ n t 10 (m/s)' 4 ), and u^ is a scale phase velocity (^6 m/s). 


*e> "-T2 

9irf ((0/f) + 1) 1 

2tT CO ^ 

describes the angular variation normalized so that f j F(e) sine d0 d<b = 1, 


o o 

where 9 is the zenith angle of k, and f = f/n, where f is the inertial 
frequency. F( e) / (4/ 9 irf 2 ) = ( e/f ) 2 /(( e/f ) 2 + l) 7 ^ is plotted in the figure. 


The maximum lies at 9 (°) = (*5/2)f = ( v3/2)(sin(latitude)/120n(rad/s)). 
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Figure 1. The angular dependence of a(k). 
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For n = 1 x 10 (rad/s) (troposphere) at Arecibo (18.3°), Sunset (40°), and 
Poker Flat (65°), 0 max = 0.23°, 0.46°, and 0.50°, respectively. 

With this form of E^(}£.) 

M 2 Qtt -9 1 / 2 E * 

® ~3T2 

u * 


and the power reflectivity is 

| p | 2 “ 4 // o(k)) g 2 ( 0 1 ,4> * ) sin8 d8 d<h 

where g(8 , ,(j) , ) is the normalized antenna gain function. If g is broad compared 
with F, as it usually is, then near the zenith the integral = 1. 

This model can be tested by comparison with measurements of ®(^) or |pP 
as functions of the parameters. Two reservations must be kept in mind, however. 
First, it has been assumed implicitly that E* is a universal constant, 
independent of latitude, altitude, etc. In fact, E* could depend on latitude 
through f and on altitude through the atmospheric density p, as it should 
theoretically. Observations so far have been inadequate to describe such 
variations. Nevertheless, in order to avoid this uncertainty, tests of the 
model should be made at a given latitude and over a small altitude range. 

Second, although the buoyancy wave model fits observed spectra rather well 
down to vertical scales as small as 9,20 m, the limit of the observed spectra, 
at some scale not much smaller than 20 m the spectrum should start to become 
increasingly attenuated by K-H instability and eventually by viscosity. Thus, 
the model can be quantitatively correct only for radar frequencies smaller than 
about 7.5 MHz (corresponding to 20 m). Unfortunately, the only suitable 
observations of Fresnel reflectivity are in the attenuated range at frequencies 
between 40 and 54 MHz (^ ra da r /2 between 3.7 and 2.8 m). In this range the 
observed reflectivity should be much smaller than the model reflectivity and the 
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dependence of k should be much stronger than k~l/2. The dependence on n, f, 
and 6 may also differ from the model. 

Nevertheless, comparisons at frequencies between 40 and 54 MHz are of 
interest. Results so far are: (1) The width of the model angular spectrum 

convolved with the radar beam is roughly consistent with the observed widths 
(GREEN et al., 1981; ROTTGER et al., 1981). (2) The model angular width varies 

inversely with the Brunt-Vaisala frequency, so that the width in the strato- 
sphere should be about 1/2 the width in the troposphere, consistent with 
observations (ROTTGER et al., 1981). (3) The model dependence on n is given by 

Mr = n^’ , roughly consistent with observations (GAGE et al., 1981). 

(4) The model ]p|^ for X ra£ jar = 7.4 m is much larger (by a factor of 100) than 
the reflectivity reported by GAGE et al. (1981), as is to be expected. 

These comparisons are satisfactory at the present level of development of 
the model and the experiments, but further, more critical, tests are clearly 
needed. Observations with radar half-wavelengths > 20 m (f < 7.5 MHz) are 
clearly needed. The angular variation might be measured by means of an inter- 
ferometer. In spite of the limitations of the model for half-wavelengths in the 
lower VHF range, further comparisons should also be made there, since the 
attenuation of the spectrum may not have first-order effects on the dependence 
on n, f, and 8. 

It should be noted that if the model is substantiated, then radar 
observations in the attenuated range yield information about the spectrum of 
buoyancy waves at small scales that is very difficult to obtain by other means. 
The buoyancy wave spectrum in this regime has implications on the cascade of 
energy in the buoyancy wave field, on the turbulent energy dissipation rate, and 
on turbulent mixing by the K-H breakdown of buoyancy waves. 

All of the proposed mechanisms for the generation of the Fresnel 
refractivity irregularities, including the present one, depend upon speculative 
assumptions. The present model is distinguished from the others by being 
quantitative and therefore testable, at least under some conditions. 
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